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Abstract. The general picture of how thermal AGNs work has become clearer 
in recent years but major observational puzzles threaten to undermine this 
picture. These puzzles include AGNs with extremely asymmetric emission line 
\^ . profiles, inconsistent multi- wavelength variability, rapid apparent changes in the 

sizes of emitting regions and in the direction of gas flow, a curious insensitivity 
of gas in some narrow velocity ranges to changes in the ionizing continuum, and 
differing dependencies of polarization on gas velocity. It is proposed that all 
these puzzles can readily be explained by off-axis variability. 
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1. INTRODUCTION 



AGNs can be divided into "thermal AGNs" and "non-thermal AGNs" (see 
Antonucci 2011) depending on whether the dominant continuum emission is ther- 
mal or non-thermal. Thermal AGNs have high accretion rates while non-thermal 
AGNs have accretion rates many orders of magnitude below the Eddington limit. 
I will argue here that we now have a reasonably secure picture of what the inner 
regions of a thermal AGN are like, but that there are many hard-to-explain ob- 
servations. I will discuss how continuum variability of thermal AGNs has to be 
non-axisymmetric on short timescales and then outline how this off-axis variability 
readily explains the hard-to-explain observations. 



2. THE STANDARD MODEL OF THERMAL AGNS 



Our current paradigm of thermal AGNs (and probably non-thermal AGNs as 



well) is that the energy ultimately comes from accretion of matter onto super- 
massive black holes. The accreting material forms an accretion disc (Lynden-Bcll 
1969). As is reviewed in Gaskell (2008), the emission of thermal AGNs is domi- 
nated energetically by thermal emission of the so-called "big blue bump" (BBB). 
A small amount of the BBB radiation is Compton up-scattered to give hard X-ray 
emission, and some of it is absorbed by dust grains and re-radiated in the infra-red. 
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Fig. 1. A schematic view of the BLR and torus of an AGN in a plane through the axis. The 
figure is approximately to scale (except that the black hole is shown too large!). Figure from 
Gaskell, Goosmann, & Klimck 2009. 



The BBB radiation is also reprocessed in the broad emission lines that are almost 
certainly present in all thermal AGNs. 

I believe that you do not understanding something unless you know what it 
looks like. Although it was (or should have been) clear from Curtis (1918) and Jen- 
nison & Das Gupta (1953) that AGNs have axial symmetry, models of AGNs were 
mostly spherically symmetric until the late 1970s. A paradigm shift took place 
between 1978 and 1980 when Blandford & Rees (1978) proposed that BL Lac ob- 
jects were AGNs seen with their beams aimed at us, and Keel (1980) discovered 
that Seyfert 1 galaxies (lower luminosity AGNs where we can see the broad-line 
region, BLR) are preferentially face on. These discoveries led to what are called 
"unified theories" of AGNs. In these theories AGNs are fundamentally similar but 
look different depending on their orientation because of either relativistic beaming 
( "beaming unification" ) and/or obscuration by a dusty torus ( "obscuration unifi- 
cation"). These theories are thoroughly reviewed by Antonucci (1993, 2011). It 
has long been recognized from the temperature of the dust that the radius of the 
torus is somewhat greater than the radius of the BLR. This is confirmed by re- 
verberation mapping measurements (Clavel, Wamstecker, & Glass 1989; for more 
recent results see Suganuma et al. 2006 and Gaskell et al. 2007). Gaskell, Klimek, 
& Nazarova (2007; GKN) argue that the BLR is simply an inner extension of the 
dusty torus and show that the BLR is flattened, highly stratified in ionization, 
and must have a similar covering factor to the torus. For more details see Gaskell 
(2009). A cartoon of the BLR and torus is shown in Fig. 1. 

The kinematics have to be consistent with this structure. As proposed by 
Osterbrock (1978), the dominant motion in AGNs is rotation and then there is a 
substantial vertical component of velocity (what Osterbrock called "turbulence"). 
Gaskell & Goosmann (2008) argue that there is also a significant amount of inflow. 
This geometry and kinematics not only explains double-peaked "disc-like" broad- 
line profiles, but also, contrary to what is widely believed, the model also readily 
explains centrally-peaked "logarithmic" profiles (see Gaskell 2010). A relatively 
small change of viewing angle is all that is needed. This is illustrated in Fig. 2. 



OFF- AXIS VARIABILITY OF ACTIVE GALACTIC NUCLEI 



279 




-6000 



-4000 



-2000 2000 

Radial Velocity (km s 1 ) 



4000 



6000 



Fig. 2. The dot-dashed line shows the Balmer line profile that would arise from the GKN 
model if it were centrally illuminated and viewed from 30° off axis. The profile with short dashes 
is the same model viewed from face-on. The thin line superimposed on this is a logarithmic 
profile. Figure from Gaskcll (2010). 



3. CHALLENGES FOR THE CURRENT PARADIGM 

While the GKN model is consistent with unified models of AGNs and explains 
broad-line strengths and profiles, there are many puzzling observations which have 
no ready explanation in the current axisymmetric paradigm (Gaskell 2010). These 
include the following: 



1 . A large fraction of AGNs have asymmetric broad line profiles and some have 
extremely asymmetric profiles with the emission line peaks displaced from 
the rest of the host galaxy by thousands of km s _1 (Gaskcll 1983; Gaskcll 
1996; Eracleous & Halpern 2003). 

2. Multi- wavelength monitoring shows that different spectral regions often vary 
independently and in a way that is inconsistent with reprocessing of higher- 
energy radiation (see Gaskcll 2006, 2008 and references therein). 

3. In reverberation mapping the discrepancies between the line responses and 
continuum variability exceed the observational errors and different contin- 
uum variability events give different time delays (Maoz 1994). 

4. Velocity-resolved reverberation mapping can show discrepant kinematic sig- 
natures (Dcnney et al. 2009) and there can be rapid apparent changes in the 
apparent direction of gas motion on timescales that are much too short to 
be real (Kollatschny & Dietrich 1996). 

5. Root-mean-square spectra can show that the most variable part of a line 
profile can have a surprisingly narrow range of velocity (see Fig. If of Denncy 
et al. 2010, for example). 

6. Broad line gas in narrow velocity ranges frequently shows puzzling insensitiv- 
ities to changes in the ionizing continuum (Sergeev et al. 2001). This is shown 
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Fig. 3. The correlation of broad-line flux with continuum flux as a function of velocity for 
NGC 4151 in six separate observing seasons. The solid line shows correlation coefficients for Ha, 
the dashed line for H/3, and the dots for Hell A4686. Adapted from Sergeev et al. (2001). 

in Fig. 3. It can be seen that each observing season there are narrow veloc- 
ity regions showing significantly weaker correlations or even anti-correlations 
with the continuum flux. 

7. The polarization changes systematically across broad lines, but the changes 
imply different kinematics in different objects (see Axon et al. 2008). 



4. OFF-AXIS VARIABILITY 

Contrary to the assumptions of the axisymmetric paradigm, continuum vari- 
ability has to be non-axisymmetric (Gaskell 2008). The reason for this is easy to 
understand. First let us consider an axisymmetric accretion flow. The idealized 
radial temperature gradient produced in an accretion flow when the only heating 
comes from the potential energy liberated at a given radius has long been under- 
stood and the resultant spectrum calculated (Pringle & Rees 1972; Shakura & 
Sunyaev 1973; PRSS). In reality there will be additional heating at a given radius 
because of the flow of heat from the hotter inner regions. This flattens the tem- 
perature gradient. The left panel of Fig. 4 illustrates the temperature gradients 
resulting from pure local energy generation (the case considered by PRSS) and 
central illumination only; the left panel shows the resulting spectra. Note the 
strong sensitivity of the spectral shape to the temperature gradient. This means 
that the temperature gradient is well constrained by the observed spectral shape. 

The spectra in Fig. 4 were calculated by assuming that the gas at each radius 
radiates as a black body and adding up the relative contributions. This is illus- 
trated in left panel of Fig. 5. For purposes of illustration the smooth temperature 
gradient of the dashed curve in the left panel of Fig. 4 has been approximated 
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Fig. 4. (a) The temperature gradients in an accretion flow corresponding to purely local 
energy generation (the case considered by Pringle & Rees 1972 and Shakura & Sunyaev 1973) 
and to central irradiation only. The dashed line shows the temperature structure deduced from 
the observed spectral energy distribution of an AGN (Gaskell 2008 for details), (b) The spectra 
produced by truncated discs with the theoretical temperature structures shown in the left panel. 
The dotted line shows the classic F v oc 1/+ 1 ' 3 accretion disc spectrum of PRSS for reference and 
the dashed line shows a F v oc v~ 1 power-law. 

as a scries of steps in temperature. The inner radius has been taken to give a 
cutoff in the spectral energy distribution around 4 Rydbergs (Mathews & Ferland 
1987). Each curve corresponds to approximately a factor of two in radius. The 
resultant integrated spectrum of the accretion flow is shown by the dashed line. 
Notice that, at a given wavelength, the flux predominantly arises from a given 
radius. For example, if we consider the flux at A3000 (~ 10 15 Hz) it can be seen 
that there are negligible contributions from the radius producing the extreme UV 
or the soft X-rays and from the radius producing the near IR. 

It has long been known that AGNs are highly variable in the optical region 
of the spectrum and at higher energies (see, for example, Fig. 4 of Gaskell 2008). 
Because the thermal emission at a particular wavelength comes from a limited 
range of radius, the variability at that wavelength must also be arising within a 
limited radius. Fig. 5 illustrates how a flare at some radius only influences a lim- 
ited spectral region. We thus expect that the further apart spectral regions are 
in wavelength, the more independent is their variability. This is indeed what is 
observed (see Gaskell 2008). The independence of variability of different spectra 
regions (see Gaskell 2008) is important because it rules out variability at lower 
energies simply being the result of reprocessing of variability at higher energies. 
Instead, there are clear cases of variability arising at lower energies without as- 
sociated high energy variability (see, for example, Fig. 7 of Gaskell 2008). Such 
variability must arise in an annulus, such as the central annulus of the cartoon 
shown in the right panel of Fig. 5. The important point is that when variability is 
limited to an annulus, the variability cannot be simultaneously around the annu- 
lus. There simply is not enough time for the annulus to coordinate its variability! 
Instead, as proposed in Gaskell (2008), the variability must be non-axisymmetric. 



5. CONSEQUENCES OF OFF-AXIS VARIABILITY 

Off-axis variability offers immediate solutions to all the problems enumerated 
in section 3. Space here only permits a listing of the solutions with brief comments. 
For more details the reader is referred to Gaskell (2010). 
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Fig. 5. (a) The theoretical spectrum arising from an accretion flow. The individual Planck 
curves (see text) are shown to scale. The dashed line is the integrated spectrum from a uniform 
disc. The thick line shows the spectrum when the temperature of the zone producing the flux 
around 1 Rydberg (the thicker Planck curve) increases by 20%. (b) Cartoon of variability within 
an annulus. 

1. Off-axis emission naturally explains displaced broad-line peaks. As an il- 
lustration of this, Fig. 6 shows how the most extreme asymmetric profile 
currently known (SDSS J0946+0139 - Boroson & Lauer 2010) can be read- 
ily fit with off-axis emission. The active region responsible for the displaced 
red peak is on the receding side of the disc near the inner edge of the region 
producing H/3. The structure and kinematics are similar to the model of 
Fig. 2. If the illumination is symmetric, the BLR shows the same profiles as 
in Fig. 2. The off-axis model also easily explains changes in line profiles as 
the active areas orbit and turn off and on (see Gaskell 2010). 

2. As has already been discussed in secion 5, off- axis variability explains why 
different spectral regions often vary independently. 

3. Different continuum variability events giving different time delays are ex- 
plained by differing locations of continuum events. Those on the far side of 
the black hole will produce shorter delays, and those on the near side will 
produce longer delays. This is quantitatively modeled in Gaskell (2010). 

4. Off-axis variability explains how different kinematic signatures can be seen 
in velocity-resolved reverberation mapping and how different events can give 
different kinematic signatures in the same object. A more rapid response of 
the blue side of a line simply means that the region of continuum variability 
is on the approaching side of the disc, rather than that the BLR is outflowing. 



5. Off- axis variability produces strong variations in the parts of the disc closest 
to it. Because different regions of the disc produce different radial veloci- 
ties, this explains how the most variable part of a line profile can have a 
surprisingly narrow range of velocity. 
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Fig. 6. The observed SDSS H/3 profile of SDSS J0946+0139 (dots) compared with the off- 
axis illumination model (thick solid line) in which the GKN BLR has been illuminated from an 
active region located at the inner edge of the H/3-emitting region. The NLR contribution to H/3 
to the observed spectum has not been subtracted. The narrow emission line at +6000 km s _1 
is [O III] A 4959. The dot-dashed line at the bottom is a similar off-axis profile to that shown in 
Fig. 1 (i.e., with central illumination and viewed from 30° off axis). Figure from Gaskcll (2010) 
where more details of the modeling can be found. 

6. Similarly, localized variability can be independent of the total flux variabil- 
ity. This explains the weak correlations or anti-correlations seen in narrow 
velocity regions in Fig. 3. 

7. Finally, off-axis emission naturally produces velocity-dependent polarization. 
Modelling shows (see Gaskell 2010) that the object-to-object differences in 
the changes in polarization across emission lines can naturally be explained 
by differing positions of the most active regions. There is no need to invoke 
different kinematics or geometries in different AGNs. 



6. DISCUSSION AND CONCLUSIONS 

Although other causes have been proposed for many of the phenomena enumer- 
ated in section 3, off-axis variability has the advantage of Occam's razor: one does 
not need to postulate different conditions in different objects. Instead, off- axis 
variability offers a simple common explanation of the phenomena. It also provides 
support for our basic picture of AGNs and for the fundamental similarity of AGNs. 
It does away with the need for exotic explanations, such as supermassive binary 
black holes Gaskell (1983). 

The off-axis variability paradigm makes specific predictions about the details 
of variability that can be tested (see Gaskell 2010 for details). The strongest pre- 
diction is that, since polarization is highly sensitive to departures from symmetry 
about the line of sight, the velocity dependence of polarization will change with 
time as different regions become active. 

Obviously off-axis variability makes some things more complicated: it severely 



284 C. M. Gaskell 



limits what we can learned from reverberation mapping and it makes looking for 
close binary black holes difficult since the profiles of a single BLR are complicated. 
On the other hand, in the new paradigm, things such as changes in the velocity 
dependence of line polarization or the velocity dependence of the correlation of line 
and continuum flux variability (Fig. 3) enable us to learn about the distribution 
and frequency of occurrence of flares close to the BLR. 
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